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low  computation  for  a more  general  body  configuration;  and  (4)  modification  of 
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I.  INTRODUCTION 


Hiis  report  describes  research  work  whose  ain  was  to  extend  a 
Magnus  force  calculation  nethod  over  rlg^t  circular  cones  to  artillery 
projectiles.  Previously,  a sethod  of  calculating  the  Magnus  force  over 
a spinning  right  circular  cone  had  been  developed  in  references  1 and  2, 
and  this  nethod  was  restricted  to  lanlnar  flow.  In  order  for  the 
methods  developed  to  be  extended  to  artillery  projectiles,  some  signifi- 
cant changes  had  to  be  made.  A listing  of  these  changes  will  now  be 
given. 

(1)  Axial  pressure  gradients  were  Introduced  into  the  boundary 
layer  equations  and  calculation  methods  developed  to  take  account  of 
the  curvature  of  the  surface  of  an  artillery  type  projectile. 

(2)  Turbulent  transport  properties  were  introduced  into  the 
boundary  layer  equations,  since  most  realistic  applications  Involve 
turbulent  boundary  layer  flows.  With  the  use  of  turbulent  transport, 
a turbulence  model  had  to  be  chosen,  and  also  new  methods  had  to  be 
developed  to  calculate  turbulent  flow.  The  major  new  method  Involved 
the  development  of  a coordinate  transformation  to  resolve  the  peculiar 
structure  of  turbulent  boundary  layers  with  an  efficient  nun^er  of 
numerical  node  points. 

C3)  The  calculation  of  the  Magnus  force  had  to  be  modified  to 
reflect  the  change  in  geometry  from  a cone  to  an  artillery  projectile. 
This  involved  changing  the  method  by  which  the  boundary  layer  forces 
were  integrated  and  changes  in  the  calculation  procedure  for  the  three- 
dimensional  boundary  laye’^  thickness. 

C4)  An  additional  study  was  carried  out  to  assess  the  feasibility 
of  using  boundary  region  methods  to  calculate  the  viscous  flow  near  the 
body.  This  study  did  not  go  past  the  stage  of  a derivation  of  the 
boundary  region  equations  for  the  artillery  projectile  geometry.  A 
detailed  description  of  each  of  these  modifications  will  be  given. 


2.  H.  A.  DujyeVt  "Three  IHmeneional  Flou  Studies  Over  a Spiwving  Cone 
at  Angle  of  Attaokf"  BRL  Contraot  Report  No.  137 1 U.S.  Amy 
Ballietio  Research  Laboratories ^ Aberdeen  Proving  Ground^  Marylandt 
February  1974,  AD  774795. 

2.  H.  A.  Divyer  and  B.  R.  Sanders,  ^Magnus  Faroes  on  Spinning  Super- 
sonic Cones,  Part  I:  The  Boundary  Layer, " AIAA  Journal.  Vol.  14, 

No.  4,  April  1976,  pp.  498-504. 


II.  AXIAL  PRESSURE  GRADIENTS 


For  conpressible  boundary  layer  flow  over  an  axisynmetric  body, 
such  as  an  artillery  shell,  defined  by  a radius  of  curvature  r(x),  the 
three-dimensional  boundary  layer  equations  can  be  written  as 

O r„  is.  + V ] 3p  9 / 9u 

*’["3x  (x-Momentum) 

^ 3x  3y  * r 3^  * ~ “'J  • “ 7 ^ ^ 3y  j (♦—  Momentum) 
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■me  numerical  solution  of  these  equations  can  be  considerably  improved 
in  accuracy  and  speed  by  transforming  coordinates,  me  coordinate 
transformations  are 
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In  terms  of  these  new  independent  variables  the  boundary  layer  equations 
become 
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The  primary  difference  between  these  equations  and  those  for  a right  cir- 
cular cone  are  the  additional  terms  involving  axial  pressure  gradients. 
The  axial  pressure  gradients  are  calculated  from  the  inviscid 

flow  equations,  which  have  been  solved  by  the  "shock  capturing"  finite 
difference  scheme.  A problem  which  aweared  was  that  the  inviscid  limit 
of  the  boundary  layer  equations  was  differenced  in  a slightly  different 
method  than  the  inviscid  flow  calculation.  When  this  occurs  the  boundary 
layer  equations  can  exhibit  overshoots  and  undershoots  as  the  inviscid 
flow  is  approached  from  the  wall.  However,  if  the  inviscid  pressure 
gradients  are  calculated  with  a procedure  consistent  with  the  boundary 
layer  finite  difference  scheme,  no  difficulties  arise. 


III.  TURBULENT  FLOW  AND  TRANSPORT  PROPERTIES 

In  order  to  form  the  boundary  layer  equations  for  turbulent  flow, 
the  dependent  variables  are  broken  up  into  fluctuating  and  mean  compo- 
nents, reference  3.  The  result  of  this  procedure  yields  the  following 
equations  in  the  physical  plane. 

(P  u r)  > (p  V r)  + Ip  (^  w)  = 0 
where 


and  the  quantities  with  a bar  on  them  imply  an  average. 

3.  J.  C.  Adamet  "Finite-Differenoe  Anatyeie  of  the  Three-Dimenaional 
Turbulent  Boundary  Layer  on  a Sharp  Cone  at  Angle  of  Attack  in  a 
Supereonio  FlaWt " AIAA  Fcaper  No.  72-186,  10th  Aeroapaoe  Soienoee 
Meeting,  San  Diego,  California,  January  1972. 

7 


p u|H+  ^ 3 r 3u 

3x  3y  r 3(J  r " ’ 37"  * W ^ ' P " 


3y  ■ 3y 


(x — Momentum) 


pu^*0V^*^^*P-!ULr'-  1 3 r 37  --T-r] 

3x  P ''  3y  r 3^  * r ~ ' 7 W w ^ - P VwM 


3y  r 3^  3x  r 3(>  j 3y  | jsy. 


(^ — Momentum) 
/.-  2" 


■ P 3y  - P v'w'  0*1^  I " P 5 (Energy) 


Tn  ordei'  to  reduce  the  number  of  unknowns  in  the  equations,  an  eddy  vis- 
cosity and  a turbulent  Prandtl  number  are  introduced.  These  quantities 
are  defined  by  the  following  relationships 
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If  the  coordinate  transformation  plus  these  definitions  are  introduced 
into  the  boundary  layer  equations,  they  become 
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The  one  further  step  to  be  taken  is  to  introduce  a semi -empirical  formula 
for  the  mixing  length,  I,  appearing  in  the  expression  for  the  eddy  vis- 
cosity. The  expression  used  was  taken  from  reference  3 and  is 


t = X 6 tanh  I — ^ ' 


X 5 


where  X = .09  and  6 is  the  boundary  layer  thickness  defined  by 
u 


y = 6 


at  ^ = .995  . 
1 


In  order  to  efficiently  calculate  a turbulent  boiuidary  layer  by  a 
numerical  method,  one  additional  coordinate  transformation  must  be 
applied.  This  transformation  allows  for  close  grid  spacing  near  the 
wall  and  wider  spacing  in  the  outer  regions  of  the  boundary  layer. 

The  transformation  is  described  in  reference  .4  and  is 


N //SNq  1/AN 

n^  = n^(K  ^ ° 


- 1)  (K 


- 1) 


1,  2,  3,  ...,  d 


where  n^  is  the  value  of  n at  the  boundary  layer  edge,  j is  the  n finite 


difference  index,  N = (j  - 1)an  and  N,  = 1.  The  values  of  K and  AN 
J J o 
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F.  G,  Blottnev,  **Vapiable  Gvid  Scheme  Applied  to  Turbulent  Boundopy 
Layersy  » Journal  of  Computer  Methods  in  Applied  Meahaniaa  and 
Engineering^  197 S.  — 
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ii 


r 


are  two  parameters  which  are  chosen  to  give  the  desired  grid  spacing. 
Acceptable  values  for  the  turbulent  boundary  layers  occurring  over  artil- 
lery projectiles  are  in  the  range  K a 1.5  and  AN^  - .05.  The  changes 

required  in  the  boundary  layer  equations  are  that  all  derivatives  with 
respect  to  n become 

\ _ i_  M 

3n  " 3N  3n 

where  3N/3n  is  evaluated  using  finite  differences  from  the  defining 
relationship  between  n and  N given  above. 

The  finite  difference  methods  used  to  integrate  the  turbulent 
boundary  layer  equations  was  the  same  as  reference  2.  A modification 
developed  not  previously  discussed  was  the  use  of  a daiiq>ing  function  to 
make  the  transition  from  the  fully  turbulent  flow  to  the  laminar  sub- 
layer near  the  wall.  The  dashing  will  be  denoted  by  D and  is  given  by 

D = 1 - exp(-y*/A*) 

I I where  the  notation  is  described  in  reference  2.  To  apply  this  correction, 

I the  mixing  length  t is  multiplied  by  D.  One  further  semi -empirical  modi- 

f fication  for  turbulence  was  also  introduced  and  consisted  of  a function 

p to  make  the  transition  from  laminar  to  turbulent  flow.  The  one  chosen 

r was  simply  a linear  increase  of  e from  0 to  its  full  value  over  four  grid 

I points.  This  modification  worked  very  satisfactorily  in  the  absence  of 

appropriate  data  on  transition. 


IV.  CALCULATION  OF  THE  THREE-DIMENSIONAL  DISPLACEMENT  THICKNESS 

In  order  to  calculate  the  three-dimensional  displacement  thickness 
A for  flow  over  an  artillery  projectile,  a new  equation  had  to  be  derived 
This  new  equation  reflects  the  differences  between  right  circular  cones 
and  artillery  projectiles,  and  is  different  than  the  one  used  by  refer- 
ence 5.  The  new  equation  is 


"I'^i’^  If " pi“i  I?  = k ipi‘^i’^1  * It  tPi^iii  ^ fc  (pi'^i’^y 


* 3t  ^^l''l®t^  * 


5.  B,  R.  Sar^e,  yheee-Dimenaioru^^  Steady,  Inviedd  Ftou  Field 
Caloulatiana  with  Application  to  the  Magmie  Problem,"  Ph.D,  Thesie, 
Department  of  Mechanical  Engineering,  University  of  Califomia- 
Davie,  California,  1974. 
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where  6^  and  6^  are  the  primary  flow  and  cross  flow  displacement  thick- 
nesses, respectively. 

It  was  also  found  that  the  method  of  solution  used  in  reference  5 
had  an  error  in  it,  and  a new  method  had  to  be  developed.  The  new  method 
uses  backward  differences  for  3A/3x  and  and  obeys  the  single  mathe- 

matical characteristic.  This  characteristic  is  defined  by  the  equation 


The  new  method  has  been  found  satisfactory  for  the  small  angle  of  attack 
cases  considered  up  until  the  present  time. 

V.  BOUNDARY  REGION  EQUATIONS 

At  large  angle  of  attack  or  whenever  the  cross  flow  velocity 
"separates, *•  the  boundary  layer  equations  as  described  above  are  not  an 
adequate  description  of  the  viscous  flow  near  the  wall.  In  order  to  des- 
cribe this  flow,  the  boundary  region  equations  have  been  proposed,  refer- 
ence 6.  The  boundary  region  equations  for  an  artillery  projectile  are 
listed  below.  These  equations  contain  additional  terms  from  the 
complete  Navier-Stokes  equations,  which  represent  viscous  effects  in  the 
cross  flow  direction.  The  equations  are  as  follows 


P i u ^ + V ^ ^ 

I 3x  3y  r 3(|i  r 
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3x  3y  " 3y 


(continuity) 
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T,  C.  Li/n  and  S,  G.  Rubirif  "A  Tujo-Layer  Model  for  Coupled  Three- 
Dimeneional  Viaaoue  and  Inoiaoid  Flow  Caloulationet"  AIAA  Paper  No. 
75- 853 t 8th  Fluid  and  Plaama  Dynamioa  Conference ^ Hartford^ 
Conneotioutt  June  1975. 
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The  solution  of  these  equations  is  beyond  the  scope  of  the  present 
task,  but  their  solution  will  be  necessary  to  fully  understand  the  Magnus 
problem  on  artillery  projectiles. 


VI.  RESULTS  AND  SUMMARY 

It  can  be  said  that  this  research  effort  has  been  very  successful. 
The  major  modifications  of  the  cone  results  such  as  axial  pressure 
gradient,  turbulent  flow  and  three-dimensional  displacement  thickness 
have  been  completed.  Dr.  Walter  Sturek  and  Dr.  Robert  Reklis  of  BRL 
are  using  computer  programs  written  by  the  present  author  and  these 
computer  programs  are  calculating  the  Magnus  force  over  artillery- type 
projectiles  for  realistic  conditions.  Dr.  Sturek  is  comparing  these 
results  with  preliminary  experimental  data,  and  the  agreement  looks 
very  good. 
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LIST  OF  SYMBOLS 


specific  heat  at  constant  pressure 


static  enthalpy 


molecular  conductivity 


mixing  length 


transformed  y coordinate 


molecular  Prandtl  number 


local  radius  of  model 


velocity  component  in  x direction 


velocity  component  in  y direction 


velocity  component  in  direction 


surface  coordinate  in  longitudinal  direction 
coordinate  perpendicular  to  local  surface 


turbulent  eddy  viscosity 


boundary  Ifyer  thickness 


boundary  layer  displacement  thickness 


three  dimensional  boundary  layer  displacement  thickness 


transformed  x coordinate 


density 


coordinate  in  circumferential  (azimuthal)  direction 


LIST  OF  SYMBOLS  (Continued) 


Subscripts 

1 edge  of  boundary  layer 

” free  stream  reference  condition 

Superscripts 

' fluctuating  quantity 

time  averaged  quantity 
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